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We derive a general formula for the paramagnetic NMR nuclear shielding tensor of 
an open-shell molecule in a pure spin state, subject to a zero-field splitting (ZFS). 
Our findings are in contradiction with a previous proposal. We present a simple 
application of the newly derived formula to the case of a triplet ground state split by 
an easy-plane ZFS spin Hamiltonian. When kT is much smaller than the ZFS gap, 
thus a single non-degenerate level is thermally populated, our approach correctly 
predicts a temperature-independent paramagnetic shift, while the previous theory 
leads to a Curie temperature dependence. 
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The nuclear magnetic resonance (NMR) chemical shift in molecular species with an open- 
shell electronic structure is mainly governed by a term known as the 'paramagnetic shift', 
a temperature-dependent term arising from the internal magnetic fields generated by the 
unpaired spin and unquenched orbital moments of the thermally populated Zeeman-split 
electronic degenerate ground state. Despite the increasingly central role played by param- 
agnetic NMR in the elucidation of the structure of metallo-proteins 1 , and in the investigation 
of the spin dynamics in novel magnetic materials 2 , only quite recently rigorous theories have 
been developed for the ab initio calculation of the paramagnetic NMR chemical shift^"-. 

Of particular relevance in this respect is the work of Moon and Patchkovskii 3 , that of 
Pennanen and Vaara 4 , and that of Van den Heuvel and Soncini^. Moon and Patchkovskii 
derived an expression for the paramagnetic shielding tensor of a spin doublet state in terms 
of its g- and A-tensors 3 . This treatment was extended by Pennanen and Vaara to arbitrary 
spin states, in the limit of weak spin-orbit coupling 4 , and later generalised by us to a theory 
that is valid for arbitrary strength of spin-orbit coupling, and arbitrary size of the degenerate 
manifold 5 ^. In Ref. |4| the authors also proposed a general formula for the paramagnetic 
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shielding tensor of a spin state subject to zero-field splitting (see Eq. (10) in Ref. 
communication we present an alternative derivation of this formula, based on the general 
theory of NMR chemical shift we have recently developed^. Interestingly, we obtain a result 
that differs from that proposed in Ref. U. In the last part of this communication we point 
out the difference, and argue for the correctness of our proposal by way of a simple example. 

We consider a molecule in the frozen nuclei approximation (also known as the 'solid 
state limit' of NMR->&). Assuming that the zero-field splitting in the degenerate ground 
state is much smaller than the energy separating the ground state from excited states and 
assuming that these excited states are not thermally accessible, the shielding tensor cr can 
be divided in two parts: cr = cr p + <x r , the first part representing the 'paramagnetic shift', 
which is due entirely to the (quasi-)degenerate ground state and which can be calculated 
from knowledge of the ground-state wave functions only; the second part representing the 
'Ramsey term', which is the only term in case of a non-degenerate ground state. These 
terms are distinguished by the fact that in the limit of vanishing zero-field splitting of the 
ground state cr r is temperature independent while er p is proportional to 1/T. The present 
paper will consider the effect of a zero-field splitting on cr p . 

The electronic Hamiltonian consists of two parts: H = H + V. Here V is the perturbation 



due to the applied magnetic field and the magnetic field arising from the magnetic moments 
of the nuclei. Ho is the Hamiltonian in the absence of these fields, but including all other 
relevant electronic interactions. This means that H includes those interactions that are 
responsible for the zero-field splitting of the ground state. If the ground manifold consists 
of u states \ip\a), eigenstates of H with energies E\, an expression of H valid within this 
manifold is: 

U) 

H = Y,Ex\4>Xa)(^Xa\ 
X.a 

Here the index A counts the energy levels of the manifold, and the index a labels a basis in 
case E\ is degenerate. For our present purpose, the term V consists of two perturbations that 
combined give rise to paramagnetic shielding: the electronic Zeeman interaction V z = — m-B, 
and the hyperfine coupling Vm = J~ ■ /*, where m is the electronic magnetic moment, B is 
the applied field, 7- is the hyperfine field and /x is the nuclear magnetic moment^. 

A general formula for the shielding tensor was proposed in Ref. la, to which the reader is 
referred for more details: 
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Here F is the electronic Helmholtz free energy of the full Hamiltonian H = H + V. Evalu- 
ation of Eq. (JI1), and retention of the temperature-dependent paramagnetic part only, leads 
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af,: = ( e^rme-^Tjdw) , (2) 



where (-)q stands for the thermal average in the canonical ensemble of H , and /3 = 1/kT. 
The expression Eq. (JSJ) can now be easily integrated 6 leading to an exact sum-over-states 
formula: 
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Here Qo — Y2x a e ~ is the partition function. 

Our aim here is to rewrite Eq. fl3]) for the case of a pure spin ground multiplet split as a 
result of spin-orbit coupling. In molecules composed of light atoms the spin-orbit coupling 
is often a small perturbation, whose effect can be treated to good accuracy in the lowest 



order of perturbation theory. It is then a well known result^ that such treatment leads to 
a spin Hamiltonian 

H = S • D • S + fi B B ■ g • S + S • A I. 
In the notation introduced in the previous section, we thus have 

H = S • D • S, mi = -/i B 9ijSj, Ti = AjiSj, 

where the indices label the Cartesian directions x, y, z and we have used /i. = giji^I to 
convert between the nuclear magnetic moment and the nuclear spin. 

The states to be used in Eq. (J3J) can then be easily found by diagonalizing H in the 
space of 25 + 1 spin basis states \S M). Naturally the eigenf unctions will depend on the 
zero-field splitting tensor D. For the moment we leave the latter unspecified and denote the 
eigenfunctions generically by \S Xa). Now applying Eq. ([3]) gives 

°fj = — — TT E 9ikAij J2 \p y2(S Xa\S k \S Xa')(S Xa'\Si\S Xa) 
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Eq. (jlj) represents the main result of this communication. 

Eq. (jlj) has to be compared with the formula proposed by Pennanen and Vaara in Ref. 
4j. Their Eq. (10) reads 

a ^ = -^TrVrY^^ijiSkS^o, (5) 

which can be written more explicitly by performing the thermal average over the eigenfunc- 
tions of H : 
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Clearly, Eqs. (jlj) and (jHj) are not equal. In fact, only in the special case D = 0, i.e. in the 
absence of zero-field splitting, the two expressions Eqs. (jlj) and (jfJJ) lead to the same formula 
for the paramagnetic shielding tensor: 
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In every other case we argue that the correct formula is given by Eq. (j3J). Note that afj can 
be expressed in a form that is only similar to Eq. (JSJ), if we take the thermal average of a 
different operator: 
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Finally, to illustrate the difference between the results presented in this communication 
and previous works, and to argue for the correctness of our proposal, we consider the simple 
but very common case of a triplet state with an axial zero-field splitting H Q = DS^, and 
axial g- and A-tensors: 
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The eigenstates of H are simply the \S M) (with S = 1), which we further denote by their 
M value alone. Thus we have two energy levels: |0) at energy 0, and |±1) at energy D. On 
evaluating the newly proposed Eq. (@| we find: 

Mb 2g ± A ± 1 - e~f^_ 

(8) 
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Previously proposed Eq. ([5]) on the other hand predicts: 

Ate 9±A± 1 + e~ /3D 
g^n kT l + 2e~P D 
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which disagrees with Eq. (jSJ) on the value of a\. Note that erf? is the same in both theories only 
because of the specific axial symmetry of this system, implying that the ZFS Hamiltonian 
commutes with the component of the spin operator along the axial direction. That the 
formula for a\ in Eq. ([9]) must be wrong can be deduced by considering the low-temperature 
limit kT <C D for D > (easy-plane ZFS anisotropy). In this situation the ground state 
is |0) and is the only populated state of the system. Therefore the shielding should be 
temperature-independent. Eq. (Q however, predicts a Curie behaviour in this limit: 
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The correct limit is obtained from Eq. (jSJ) and is indeed a constant: 

fi B 2g ± A ± 
en — >■ — . 
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